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Математичне моделювання в1дносин елшав 
у задачах оптимально] кластеризаци об’ект!в 


В статье рассматриваются конструктивные средства математического и компьютерного моделирования 
отношений (включения, пересечения, касания, непересечения) эллиптических объектов. Определяется 
полный класс свободных от радикалов аппроксимаций рШ-функций для эллипсов и их дополнений. 
Строится математическая модель задачи оптимальной кластеризации эллиптических объектов в виде 
последовательности задач нелинейной оптимизации. Приводятся результаты тестовых примеров. 
Ключевые слова: математическое моделирование, эллипсы, р№1-функция, 

оптимальная кластеризация. 


У статт! розглянут! конструктивн! засоби математичного та комп’ютерного моделювання в1дносин 
(включення, перетину, дотику, неперетину) елштичних об’екив. Визначено повний клас апроксимащй ры- 
функщй (що в1льн! в1д радикал!в) для елшсв та 1х доповнень. Побудовано математичну модель задач 
оптимально! кластеризацй елштичних об’ект!в у вигляд! посл1довност! задач нелйно! оптим1заци. 
Наведено результати тестових прикладав. 

Ключов! слова: математичне моделювання, елшси, р-функщя, оптимальна кластеризашя. 


Тье агафе сопз1Аег$ сопзгасйуе 1001$ оР таетайса! то4дейп» ап сотршег зпиШаНоп {есбтаие оЁ 
т@айопз$ (шсГазюп, поп-пиегзесипе, юисЫшт, пщегзесип®) Бебуееп еШрис оБ]ес5. А сошр!ее с1аз$ оЁ 
Нее гаФса! арргохппаНопз$ оЁ ры-Р№псйопз юг еШрзез ап4 Фет сотр!етеп аге 4еЙпеа. У\е ргом4е а 
тафетайса! по4де| оЁ фе орНта| стает» о еШрис обес ш Фе Ююгт оЁ сопзгате орипизайоп 
ргоЫет. А патбег оЁ сотршаНопа]| гези $ аге г1уеп. 

Кеу мог4$: тафетайса! тодейпе, ерзез, ры-Рапсйоп, орйта] са$епипе. 


При создании систем искусственного интеллекта, предназначенных для решения 
оптимизационных задач размещения (Раскте&Си тэ) [например, 1-4], в частности 
задач кластеризации, конструктивным средством аналитического описания отношений 
(включения, пересечения, касания, непересечения) геометрических объектов является 
метод р№1-функций Стояна [5-7]. 
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В работах [6-8] приводится полный класс рш-функций для неориентированных 
базовых объектов, а также уникальный метод генерации математических моделей, 
основой которого является построение рШ-деревьев. Поскольку класс односвязных 
базовых объектов не включает объекты, обладающих формой эллипса, целью данного 
исследования является построение рШ-функций для эллипсов. 

В работе [9] построен нулевой уровень рШ-функции (годограф) для пары эллипсов 
в параметрическом виде. Вопрос построения р-функции с учетом непрерывных транс- 
ляций и вращений является по-прежнему открытым. При этом для формирования рН1- 
функций важным является использование свободных от радикалов функций, имеющих 
степень не выше второй, поскольку методы решения оптимизационных задач размеще- 
ния используют методы локальной оптимизации. 

Цель статьи — разработка конструктивных средств математического модели- 
рования (р№1-функций) для описания отношений эллиптических объектов в анали- 
тическом виде. 


Объекты и их отношения 


Пусть имеется эллипс Е, заданный большой и малой полуосями а и Б, круг С 
радиуса г и объекты вида: Е” =Ю^\ниЕ, С" = В? ШС. 
Полагаем, что начало собственной системы координат объекта АЕ {С,Е, С — Е“} 
находится в его центре симметрии. 
Вектор и = (х, у, 0, й) характеризует параметры размещения (х, у, 0) и коэффи- 
циент гомотетии й объекта А ‚где (х, у) — вектор трансляции, 69 — угол поворота, й>0. 
Рассматриваются следующие отношения между объектами А(и|) и В(и_): 


— пересечение: иА(и| ) Г ШЁВ(и> ) = ©; (1) 
— касание: А( и) б шВ(и >) = ФиНА(и,) п НВ(и.) = ©; (2) 
— непересечение: А( и) п В(и о) =; (3) 
— включение: А(и|) <В(и>) > шШМ(и) пише В*(и >) = ©, В =В?\пиВ. (4) 


Ры-функция Ф^В [5-7] определена для пары рЫ-объектов А(и|) и В(и> ), непре- 


рывна и удовлетворяет следующим свойствам: Ф АВ > (), если выполняется условие (3), 
ФЛ = 0, если выполняется условие (2), ФЛВ <0в случае (1). 

% 
Отношение включения (4) описывается Ф В 0 
В работе [3] доказано, что произвольный рН1-объект, граница которого форми- 
руется дугами окружности и отрезками прямых, всегда может быть декомпозирован 
на базовые объекты. 


Тогда ры-функции Ф^В с { ФА, ФРС, ФСЕ } для рассматриваемых пар объек- 


тов А и В определяются как ФАВ 


= {Ф,,1=1,...,п}, где Ф; — ры-функции для 
базовых объектов, п=пд`Ив, пд(пв) — число базовых объектов, реализующих декомпо- 


зицию объекта А(В). 
В работах [7], [8] построен полный класс рш-функции для базовых объектов. 


Юй * 
В этой связи представим аппроксимации объектов Е и Е ввиде объединения базо- 
вых объектов. 
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* . $ АВ 
Аппроксимация объектов Е, Е . Для построения ры-функций Ф“”, А,Ве! 


'Е(С,С ОВ. Е} ‚ осуществим  -аппроксимацию объектов Е, Е” базовыми объектами. 

Погрешность & зависит от значений а и Б. Аппроксимируем границу эллипса 
Е дугами окружности. Построенный овальный объект Е является -аппроксимацией 
эллипса Е . При этом выполняется соотношение: Е < Е р 


г ша #20 
Объект Е представим в виде объединения базовых объектов вида Е = () В;, 
= 


(рис. 1), где В, ={С,О,К},а именно 

Е= С, 92930 4 ОК$, (5) 
С: - круг радиуса г с центром в точке (Д|,0), Сз — круг радиуса г с центром в 
точке (—Д1, 0), где г=а-Ду, ДА} = @-БЖа+ь+ а? +Ь?) ; )› -— круговой сегмент 
р> =С. То, здесь С. — круг радиуса В =Ь+А. с центром в точке 05 =(0,-А.), 
Д> = аа? +57); р — треугольник, ь =соп\у,У›,[}, две стороны 


которого образованы касательными к кругу С. в точках у, у›: у=(-Ь.6-Ар,а.6), 


2 
у› =(Б-0+Ау,а:5), а, Е’; О.А - круговой сегмент 


р. =С. ПТ. , образованный пересечением круга С и треугольника Т д где Сл — круг 
радиуса К=Ь+А. с центром в точке (0,Д4.); Ту — треугольник, уе =. 


= сопу{Уз,уд,[2}, две стороны которого образованы касательными к кругу С й В 
2 
. у р 
точках уз, У4: у; =(-Ь.6-Д|1,-а:6), у4 = Ф.б+Ау,-а-5), [, = (0, д; 
ао+А 2 
К5 - выпуклый многоугольник, К5 = сопу{У|, У, Уз, Уд}. Заметим, что полюса объ- 
ектов С, 2, Сз,Од,К5 совпадают с центром эллипса Е. 


* < < < 
Рассмотрим объект Е . Пусть Е -— эллипс, заданный большой и малой полуосями 
ау и В. Рассмотрим эллипс Е, с метрическими характеристиками а=а-= и 


* „* 
Ь=Ь, -—=. Полагаем 8 ‚ при котором выполняется соотношение Е СЕ;. 


Аппроксимируем объект Е, объединением базовых объектов вида (рис. 3): 
ре 9 об ы (6) 
где в — круг радиуса К с центром в точке (0,Д.), В — круг радиуса В с 
центром в точке (0,-Л45), К=Ь+Д.; С| — объект С =Р ПС! (рис. 2), где С!- 
круг радиуса г=а-АД| с центром в точке (Д|,0), В, — полуплоскость, проходящая 
через точки ъ, уд: у =(-Ь.6-Ау,а.6), уд = (Ь.б+Д|,-а:5), С› — объект 
С. =РПСо, где С› - круг радиуса г с центром в точке (-А1,0), Р - 


полуплоскость, проходящая через точки у2, уз: у› =(:6+Лу,а-6), уз =(-Ь-6-ЛА,-а-5). 
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Рисунок 1 — Аппроксимация эллипса Е базовыми объектами 
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Рисунок 2 — Объект С =Р ПС: 
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* 
Рисунок 3 — Аппроксимация объекта Е базовыми объектами 
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Ри-функция для объектов В, и В Пусть имеются эллипсы Е иЕ2. Осуще- 
ствим аппроксимацию эллипсов Е' и Е? вида (5). Тогда рш-функция для объектов 
Е иЕ определяется так: 

ФЕ" =шт{Ф,Фо,Фз,Ф4,Ф5\, (7) 

Ф = шт ФС! ‚Ф =шш Фр; ‚ Фз=шш фС3В: : 
В:ЕЗ В;ЕЗ В;еЗ 
Ф.А - ши Ф2481, Фу ша ФАЗВЕ, 
В:ЕЗ В:ЕЗ 

‚ФСС — ры-функция для двух кругов, ФСР —рь1- 
Ф2?Р 
ФАК 


где 3=|С7.02.с3.04,кз} 


Фо — ры-функция для двух 


— ры-функция для двух 


функция для Си О, — ры-функция для С и К, 
круговых сегментов, Ф* ры-функция для Д и К, 
выпуклых многоугольников. 

Ри-функция для объектов Е и С*. Пусть имеется эллипс Е и круг С. Осу- 
ществим аппроксимацию эллипса Е вида (5). Тогда рш-функция для Е и С* 
определяется так: 


1. * 1/.* 1.“ 1/-.* ]* 
ФЕС -шш/ФСС ФРС фбзС фР4С фЕ5 || (8) 


где ФСС - ры-функция для Си С*, ФС" = ры-функция для Ри С*, 


= ры-функция для Ки С*. 


. = * * 
Ри-функция для объектов Е и С. Осуществим аппроксимацию объекта Е 


вида (6). Тогда рш-функцию Ф СЕ можно представить таким образом: 
Е ФСГ фСС>’ фСби, фСб2 |. (9) 
сс” 22 ба 
где Ф — ры-функция для Си С*, Ф ^^” — рш-функция для Си С. 


— * ь 
Ри-функция для объектов Е и Е . Пусть имеются эллипс Е и объект Е. 


Осуществим аппроксимацию эллипса Е и объекта Е’ вида (5), (6). Тогда ры- 


дд* 


функцию ФЕЕ можно представить так: 
ЕС” : = 1 1 1 1 1 
где ФЕС — ры-функция вида (8), ФЕС = пи. ФСС, ФРС фС3б фР4б фК5б| здесь 


ФСС, ФР? © ФК — ры-функции для С’ и базового объекта С, С, К соответственно. 


Ры-функции (7) — ( 10) являются 5 -аппроксимацией рш-функций для пар 
объектов Е! и Е2, Еи С’ Си Е”, ЕиЕ” соответственно. 

Ры-функция ФВ 
мумов гладких функций. Поэтому для каждого рШ-неравенства Ф^В >0 можно но- 


строить рШ-дерево, концевым вершинам которого соответствуют системы р№-не- 
равенств с гладкими функциями. 


В общем случае является композицией минимумов и макси- 
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Математическая модель. Рассматриваются оптимизационные задачи кластериза- 
ции М непересекающихся объектов Т;={С,Е},1еГу в контейнер О=Е{С,Е} 


минимальных размеров. 
Размещаемые объекты Т;,1е Гу характеризуются переменными (х;, у;, 0;) ‚ т.е. 


допускаются непрерывные трансляции и повороты, при этом полагаем й; =]. 


Область размещения © характеризуется переменным параметром й, при этом 
полагаем х=0, у=0,0 =0. 


Пусть ие А° - вектор переменных метрических характеристик контейнера © 


и размещаемых объектов Т.,1 Е Гу, К” — арифметическое евклидово пространство. 
Математическая модель задачи кластеризации имеет вид 


шш Е(ы), {95 
иЕЙеВ° 
\ = {иЕК° : (и) > 0}, (м) = шт{Ф. (и), Ф/ (и), <=1,2,...,А,1=1,2,..., М}, (12) 


где Р(и)=Й, Ф_ (и) — ры-функция непересечения для пары рШ-объектов, Ф, (и)— 
ры-функция включения рШ-объекта в ©, ^=0.5.(М-ПМ. 

Из соотношения (12) следует, что Х(и)>0, когда Ф'(и)>0 для каждой пары 
объектов А и В, А, ВЕ{Т|,Т.,..., Тм}, АВ и Ф"(и)>0 для каждой пары объектов 


О” ид, АЕ(Ту,То,...,Тм}, О’=В*\иЕО, шЕО - внутренность области О. 
т 
Поскольку И = (м, ‚ где У", ={иЕ ^° :х,(и)> 0} , задача (11) - (12) сводится к 


К= 
оптимизационной задаче [6]: 


Е(и) = ши Е), Е(и,),...Е(и,)}, (13) 
Р(и,.)= пуп (м), К =1,2,..9 <’. (14) 


Для решения задачи (13) — (14) используется подход, основанный на методе 
ветвей и границ, с набором правил отсечения, учитывающих оценки решения сверху, 
симметрию дерева решений и вырожденность неполных систем. 

Численные эксперименты. В качестве тестовых примеров рассматриваются задачи 
оптимальной кластеризации объектов в круге и эллипсе (рис. 4). 


2=С 


(а) (6) 


Рисунок 4 — Размещение объектов в области, соответствующее точке локального 
минимума и =(х,,у,,6,,..., ху, У, х,№ ), (а) пример - 1; (6) —- пример 2 
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Пример 1. Исходные данные: ©=С, Т, =Е|, Т> =Е2, Тз = Сз, 
и. =1, и=(ж, У1, 901, х2, у2, 02, хз, уз, 05, Й), а=4, В =3, а,=5, Б,=3, г, =3. 


Точка локального минимума: и’ = (-2.667336, -2.498968, -12.057380, 3.696938, - 
0.113279, -4.681757, -1.620644, 3.494998, 0, 6.852467).Значение функции цели в точке 


и’ = (х, а у. 9 ый.) ‚1 =Е(и’)=6.852467. 
Пример 2. Исходные данные: ©9=С, Т, = Е, Т› =Со, 
еб боев те 
Точка локального минимума: 

и’ = (-2.371988605, 1.625690207, -0.540419525, 3.482693316, -1.121795491, 0, 1.760897765). 
Значение функции цели в точке м” =(х,у,,0,,... хи, ун,0л,й”):й =Еи)=1.760897765. 


Для решения задач локальной оптимизации вида (14) используется 
ТРОРТ [10]. 


Выводы 


Приведенные рш-функции позволяют описывать математические модели задач оп- 
тимальной кластеризации эллиптических объектов в виде задач нелинейной оптимизации 
и применять для их решения известные методы математического программирования. 
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ВЕЗОМЕ 


А.У. РапКгаоъ, Т.Е. Котапота, 1.А. ЗиББоа 
Матетанса Моаейиз о} Кеаноп5 ор Е|руе5 
т Орита[ ОБлес! Сияеттз 


Тре си@те ап раскше рго ет 1$ а рай оЁ сотршанопа! геотешу Фа Ваз псь 
аррИсайоп$ ш гаптепё шдизгу, зВееё теа| сит, Випиоге таК1пе, зВое тапиасвийио, 
ес. ТЬе соттоп (азК 1$ ю р!асе а сепат зе оЁ Нгиге$ оЁ зресйе4 зВарез апа $17е5 ут а 
отуеп агеа (сотезроп4те ю а Нее оЁ теа1, а зар оЁ{ехШе, ег.). Тве обес тия Бе 
раске4 тю а <1уеп сощашег ог сиё оиё оЁа хлуеп реесе оЁ тажепа]|. То пишитияе уа%е, ог 
шшияе Фе зрасе изе4, ог тахи17е фе питбег оЁ обес, опе маш ю расКк Фе 1а(ег аз 
ОоВИу аз роз$1Ые. ис ргоЫепл$ аге Кпо\г аз си@пе ап4 раскте (С&Р) ргоЫетлз$. ш зоте 
аррИсайоп$ Феу аге саШе пезип® ргоептз, тагКег таКтх, зюск сито, сощмаштен 
ргоЫетл$, ес. Тре ргоШет 15 МР-сотр!ае, ап4 а$ а гезай пеайу а ргасйса| а]гоп11$ 
и6117е Беипт$Нс$ ап аге гезблие4 1ю обес оЁ сецаш зппр!е Варез, ап ппрозе уапои$ 
Ппкайоп$ оп Фет 1ауойё. ш Го 4итеп$1оп$, поз те@о4$ \огК ми роухоп$ оту, офег 
Варез аге зппр[у арргохипае4 Бу роугоп$ (а поаЫе ехсерйоп Беше \мсВ аПо\у$ сисаг 
зВарез). Оес($ изиаПу Вауе а Нхе опещайои, 1.е., Феу саппоЕ Бе Нееу гоаже4. ТВе агафе 
соп$1Аег$ сопзиласйуе (0015 оГ тафетайса! тодейпх ап сотрщег зпаайоп {есбтачце оЁ 
те]аноп$ (псзюоп, поп-пщегзесипе, юисШие, ищегзесипе) Бебмееп гобайпе еШрис обес. 
А сотр!ее с1аз$ оЁ га@са|-Ёгее арргохипаноп$ оЁ рШ-РапсНоп$ Гог еШрзез ап ШФег 
сотретеп(5 аге 4ейпе4. У!е ргоу14е а тафетайса| по4де оЁ Фе орита| сауение оЁ 
еШрис обес ш Фе юпт оЁ сопзгаше4 орипизайоп ргоет. А питбег оЁ сотршабнопа| 
тези!5 Гог орйта| сазепи» о еШрис обес аге слуеп. 


Статья поступила в редакцию 16.07.2012. 
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